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Scattering theory in the Gell-Mann and Goldberger formulation is slightly extended to render a 
Hamiltonian quantum mechanical description of the neutrino oscillations. 
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I. INTRODUCTION 

The concept of neutrino oscillation has become very fa- 
miliar to physicists during recent years and the subject is 
discussed in numerous articles. Nevertheless, there does 
not appear to be a Hamiltonian treatment which pro- 
ceeds from fundamental quantum-mechanical principles. 
The present work represents an attempt to fill this deficit 
starting from the Gell-Mann and Goldberger formulation 
of scattering theory [1]. There are very few new results, 
but the interpretation is somewhat new. The theory is 
limited to scattering in vacuum and does not include mat- 
ter effects [2, 3]. 

The neutrino oscillations have been identified as a po- 
tential source of information about fundamental aspects 
of particle theory by Pontecorvo [4]. Bilenky and Pon- 
tecorvo provided an initial theoretical analysis of the 
neutrino oscillation as a subject of its own merits [5]. 
Kayser [6] introduced wave packets for neutrinos. Rich [7] 
offered a space-time approach where neutrino emission, 
propagation, and absorption are treated as a single pro- 
cess. Giunti, Kim, Lee, and Lee [8] further developed an 
approach based on Feynman diagrams and wave pack- 
ets. Recently, Akhmedov and Kopp [9] summarized the 
current theory status. 

Davis, Harmer and Hoffman reported a deficit in the 
flux of neutrinos from the sun in 1968 [10]. Major collab- 
oration reports on measurements of the neutrino oscilla- 
tion are available from Super-Kamiokande [11], SNO [12], 
and KamLAND [13]. Most recent results of T2K and MI- 
NOS are described in Refs. [14] and [15], respectively. 

Section II introduces the required elements of the for- 
mal scaterring theory. The neutrino oscillations are dis- 
cussed in Section HI, where it is explained how the stan- 
dard formula emerges in their theory. Section IV con- 
cludes the paper. Appendix A explains relevant details 
of perturbative energy denominators. 



II. 



FORMAL THEORY 



The formal scattering theory of GcU-Mann and Gold- 
berger [1] starts form a full Hamiltonian H which is a 
sum of a free part Hq and the interaction part Hj. An 
incoming state \'^i{t)) of energy Ei is gradually built over 



time r = e ^ according to the formula 

\^.M^ef dTe'^e-^^(*-^)|$,(r)) (1) 

= e-^'^* -10,). (2) 

In this formula, the state |0i) is an eigenstate of Hq, 

Ham - ^.10.) , (3) 

so that 

|<I>,(T))=e-^^-^|</..). (4) 

The norms of states are not changing with time. The 
probability that the system is in a state [$/) at time t is 



UJf,{t) = 



where the amplitude Afi{t) is 



Af^{t) = (0/1 



E^- H + ie' 



(5) 



(6) 



The differential cross section for the transition i — >■ / is 
equal to the transition rate, 

PMt)^j^^Mt), (7) 

divided by the flux of incoming particles. 

Gell-Mann and Goldberger argue that for i ^ t 

\AMt)\^^\AMOW (8) 

and they focus on the case of i = 0. However, the time t 
in the neutrino oscillation experiments is typically greater 
than T. For example, in the T2K experiment, the initial 
pion states are generated from a carbon target of size of 
about 1 m (using energetic protons) and the pions only 
move about 100 m before Hj nearly certainly turns them 
into states of neutrinos and muons. This means that r 
is shorter than about ^ 1 /xs. Since the initial neutron is 
located about 300 km away, the time t must be greater 
than 300 km/c ^ 1 ms. The time t is thus expected to 
be more than 10'^ times greater than r. 

In addition, the formal scattering theory requires that 
the linear dimension of the region in which the states are 
normalized is much greater than r times the group velic- 
ity of the incident particle wave trains. This condition is 
satisfied in the T2K experiment when the required region 
contains both Tokai and Kamioka. 
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To describe scattering in these circumstances, the for- 
mal scattering formula needs to be extended to large 
times t. Using notation 

^U*^ = p F ^J ^f\ ^"''''-"^'I'Pi)' (9) 
tji — tjf + ie 

Rf,{t,e) = (^/|if,e^(^/-^)*-^— |— -|<^,) (10) 

= {^f{t)\Hj\^i{t))/ (11) 

one obtains from Eq. (6) that 

^fiit) = 5%{t) + lr—-Rfi{t, e), (12) 

Afi{t) = -iRfi{t,e). (13) 

In the spirit of [1], one can suggest that Eq. (12) is useful 
because it exhibits the energy dependence of A fi (t) in the 
vicinity of Ef = Ei. The resulting transition rate is 

\Af;{tW = 2lm[6fi{t)Rf,{t,e)] 

+ (^.-|)^ + e^ '^^-(^'-^'^- ^''^ 

Since in the neutrino oscillation experiments such as T2K 
the time r is shorter than t, one cannot send e = 
to for finite t. Instead, e approximately accounts for 
the experimental energy uncertainty and determines the 
width of density in energy of final states / at energy E^. 

The parameter e also smoothes out all contributing 
amplitudes as functions of energy. As a result of this 
smooting, the interference effects in the complete tran- 
sition rate are described by the standard oscillation for- 
mula, see Eqs. (28) and (29). In the remaining part of the 
article, it is shown how this happens when one approxi- 
mates the full Hamiltonian H in the exponential factors 
in Eqs. (9) and (10) by Hq. 

Further analysis will involve specification of Hj. We 
choose to continue focusing on the example of T2K. 

III. NEUTRINO OSCILLATION 

In the T2K experiment, neutrino oscillation occurs in 
the process 

7r"'"n — )• /i"'" /U~ p , (15) 

where the initial state consists of a 7r+-meson prepared 
in Tokai and a neutron, n, located in Kamioka. The 
final /x"*" is produced in decay of tt"*" in Tokai, and /tx~ 
and proton, p, emerge in Kamioka. The distance covered 
by the intermediate neutrinos, which are created in the 
interaction responsible for tt "'"-decay and annihilated in 
the interaction that changes n to p and creates /i^, is 
practically equal to the distance from Tokai to Kamioka. 
We assume that the only particle detected in the final 
state is Detection of other particles does not change 
our conclusions. 



Since the dominant interaction processes are tt"*" — >■ 
^'^Vjx and nv^ — > p/i^ and the neutrinos carry momenta 
and energies on the order of 1 GeV, the most appropri- 
ate interaction Hamiltonian density for obtaining Hj by 
integration over space is [16-18] 

G 

^/ = cosi?c M7"(l - 15)^^ P7a(l - 9Al5)n 

- i-^ I'ul^ii - l5)lJ- da-K^ + h.c. . (16) 
v2 

The resulting Hj is translation invariant and hence con- 
serves three-momentum. The interaction does not con- 
serve the eigenvalues of Hq, which are sums of energies 
of the form E = ^JwP- -\-p'^ for any particle of mass m 
and three-momentum p. 

The central issue is that the neutrino states that are 
created or annihilated by Hj are not eigenstates of Hq. 
One can consider Hi at the hadronic level of Eq. (16), 
or at the level of electroweak bosons that interact with 
quarks in tt"*" and nucleons. In both cases, the neu- 
trino states that are created or annihilated in the in- 
teractions with muons, called /x-neutrinos and denoted 
by f^, are thought to be properly described by a unique 
combination of (most likely) three quantum fields Vi{x), 
i = 1,2, 3, each of which corresponds to neutrino states 
that arc; ciigenstate of Hq and having different masses. In 
the case of p^, 
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fc'^(x) =^[/^ii/i(a;). (17) 

i=l 

Recent data on the mixing coefficients, U^i, can be found 
in [19]. 

The standard formula for neutrino oscillation in the 
quantum process (15) is obtained from Eq. (14) when 
the Hamiltonian H in the exponential factors in Eqs. (9) 
and (10) is approximated by Hq. In addition, one also 
expands {Ei — H + ie)~^ in Eq. (10) up to second power 
of Hj. The result is 

iiEf-Ho)t 

Rfi{t,e) = {^f\Hi- ^^Hi\cj,,) + --- , (18) 

where \4>i) = |7r+n) and \4>f) = j/Lt"*" /x~ p). 

There are only two types of intermediate states that 
contribute to Rfi between the operators Hi in Eq. (18). 
The first type contains a neutrino, /Lt+, and n. The second 
type contains an anti-neutrino, tt"*", ^,~ , and p. Only 
the first type of intermediate states can lead to a small 
energy difference in denominator, due to Ei — Hq + ie. 
Thus, the transition rate in Eq. (14) is dominated by the 
contribution from the first type of intermediate states 
and the second-type contribution can be neglected. 

Namely, the energy-denominator associated with an in- 
termediate state with a neutrino of mass rrii is 

D, = E,- E,,^ + ie , (19) 

where E,^ denotes the energy transfer from tt"*" to n, 

E, = E^+ - E^+ , (20) 
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and the neutrino energy eigenvalue of Hq is 



(21) 



All eigenvalues ascribed by Ho to all intermediate states 
with neutrinos contain E^. with one and the same phys- 
ical momentum transfer Pi, irrespective of the neutrino 
mass. The value of is implied by the three-momentum 
conservation in the translation-invariant Hamiltonian in- 
teraction terms. 

While the real part of Di in Eq. (19) must be on the 
order of e or smaller according to the density in energy 
exhibited in Eq. (14), the energy denominators associated 
with the intermediate states with anti-neutrinos have real 
parts on the order of —2E,y. Therefore, the second-type 
contributions to the transition rate are about (c/E^)"^ 
times smaller than the first-type, and are neglected. 

Consequently, Rfi in Eq. (18) is a sum of three am- 
plitudes each of which corresponds to the intermediate 
state with a virtual neutrino of a different mass. When 
one evaluates the time derivative of the modulus squared 
of the sum at time t = L/c, one obtains 9 terms which 
together combine to 



di 



t=L 



(27r)3,5(3)(p^+p„-p'^+ 
2e 



p-pV) (22) 



{E^ + E^-E^+-E^--Ep) +e2 

in a self-explanatory fashion. The factor 1^!^]^ results 
from Hj. The factor is 



where 



Gji 



1 



1 1 



AE,,E,^ D*Di 



(23) 



(24) 



The energy denominators Di are different for different 
values of i if the masses rrii in the eigenvalues i?^. are 
different. The differences between the eigenvalues are 
very small if the mass differences are very small and if 
the masses arc very small in comparison to the physical 
momentum transfer \p^\. Monte-Carlo studies of T2K 
Collaboration suggest that these conditions are satisfied. 
Under the condition that 



2|P.| 



(25) 



one obtains the approximate Gji that does not depend on 
the neutrino masses. This is shown in Eqs. (Al) and (A2) 
in Appendix A. Taking the common approximate factor 
Gji out of the sum of 9 terms, one arrives at 

= . (26) 



This result implies that the rate of counting /i~ in a dis- 
tant detector depends only on the ratio L/\piy \ as in the 
standard neutrino oscillation formula. 

Namely, the numerator in Gi, in Eq. (26) is the stan- 
dard, distance-dependent ^"-detection probability. 



1 - sin^ (26*23) sin' 



^{mym^,)L/{2\ff.\) 



2 AmjgL 



(27) 



(28) 



More precisely, the ratio of counting-rates at two dif- 
ferent distances between the neutrino detector and a 7r+- 
sourcc, such as L/ar ~ 300 km and Lnear ~ 280 m in the 
T2K experiment, is 



t=L, 



i{Lfnr) 



Pa 



t i-^near 



(29) 



This is how the perturbative scattering theory explains 
the standard neutrino oscillation formula. 



IV. CONCLUSION 

A simple extension of the Gell-Mann Goldberger for- 
mulation of scattering theory to the case of a long base- 
line experiment, leads to the standard neutrino oscilla- 
tion formula that describes the distance-dependent ra- 
tio of muon-counting rates via Eq. (29). This conclu- 
sion holds provided the condition (25) for a pion-beam 
preparation time r = is satisfied and one neglects 
all interaction effects between the initial pion decay and 
neutrino absorption in a detector. This is formally facili- 
tated by replacing H with Hq in the exponents in Eqs. (9) 
and (10). 

The oscillation formula is a result of the interference 
between scattering amplitudes mediated by virtual states 
with neutrinos of different masses. All these states medi- 
ate a transfer of the same physical three-momentum p,^ 
and energy E,y. However, their contributions to the total 
amplitude differ as functions of and E^. According 
to Eq. (23), only a sufficiently large e can smooth out 
the differences between factors Gji and yield Eq. (26), so 
that the net effect of the interference can be described 
using the standard oscillation formula. 

The condition (25) indicates that the time of prepar- 
ing the beam of tt^ mesons must be sufficiently short, 
the neutrino momentum must be sufficiently large, and 
the neutrino masses sufficiently small for the oscillation 
to occur in agreement with the standard formula. Other- 
wise, one might expect deviations, whose details can be 
deduced from Eq. (23). 

Regarding the replacement of H by Hq in the expo- 
nents in Eqs. (9) and (10), one may observe that precise 
calculations of weak interaction effects following from Hi 
of Eq. (16), require understanding of the cutoff depen- 
dence that appeared already in Pontccorvo's work [4]. 
One also ought to consider interactions with matter [2, 3]. 
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It should be stressed that the conclusion concerning the 
ratio of transition rates does not automatically translate 
to the ratios of entire cross sections. Evaluation of cross 
sections includes averaging over incoming and integra- 
tion over outgoing states and may involve wave packets, 
density matrices, entanglement, and various experimen- 
tal cuts. Such evaluation is a formidable task in a fun- 
damental theory. For example, if neutrinos are coupled 
through intermediate bosons to quarks, and one attempts 
to evaluate the complete scattering process including the 
quark structure of hadrons, in addition to the nuclear 
binding effects, an entire host of additional theoretical 
issues arise. To be more specific, the Fermi motion of 
quarks in nucleons, or of nucleons in nuclei, is a subject 
of study in its own right [20] and currently available level 
of theoretical analysis could certainly be improved, e.g., 
see Eq. (2) in [14]. 

Finally, we note that the Gell-Mann and Goldberger 
discussion of a connection between the formal scattering 
theory and S'-matrix formalism in the interaction pic- 
ture involves the limit of e — >■ 0. This limit requires 
that the beam preparation time r is much longer than 
the period in which interactions causing scattering may 
happen. If this condition is not satisifed and one can- 
not take the limit of e — >■ 0, one has to proceed from 
fundamental quantum-mechanical principles rather than 
taking advantage of simplified formulae in which e — > 0. 



Appendix A: Energy denominators 



Factor Gji in Eq. (24) can be rewritten as 



Gji 



1 



1 



1 



(Al) 



which is a valid approximation when \Ei,^ — Ei,.\/2 <C e. 
In addition, using condition (25), one obtains 



G 



^E,-\p,\y + e^ 



(A2) 



which does not depend on the neutrino masses and leads 
to Eq. (26). 

Eq. (Al) differs from Eq. (A2) by the terms that 
in a leading approximation are inversely proportional 
to the neutrino momentum and directly proportional to 
the heaviest neutrino mass squared times 2A/(A-^ + e^)j 
where A = E^ — \pv\- Therefore, the leading correction 
to the standard oscillation formula due to using Eq. (26) 
instead of Eq. (23), may be expected to be smaller than 
ct/L, which currently means smaller than lO'^ in T2K. 
In any case, its actual size depends on the pion beam 
preparation and the Monte Carlo simulations may even 
average it to nearly due to the varying sign of A in data 
sampling. A special sampling that secures only inclusion 
of cases with one sign of A would be required to facilitate 
studies of this correction. 
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